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1 . I ntrodu  ct i on 


In  the  report  (Moritz,  1976)  we  have  seen  that  the 
local  behavior  of  a covariance  function  (c.f.)  of  the  gravity 
anomaly  A g can  be  character i zed  quite  well  by  means  of 
three  constants: 


C ....  variance  of  A g , 

0 

5 ....  correlation  length  (argument  for  which  the 

c.f.  has  the  value  CU)  = j C,)  . 

G ....  variance  of  the  horizontal  gradient  of  Ag  . 


Instead  of 
related  to 


we  may  also  use  the  curvature  parameter  \ 

G by  (ibid.,  p . 25) : 

0 


X 


G E*/C 

o o 


(1-1) 


As  regards  the  global  behavior  of  a covariance 
function,  the  principal  requirement  is  to  fit  the  lower  degree 
variances  c , say  from  c,  to  c n . , as  determined  from 
satellite  and  gravimetric  data. 

As  an  example  of  a set  of  local  parameters  for  a 
global  c.f.  we  may  take  the  rounded  values 

C = 1800  mgal“  , 

o 

5=80  km,  (1-2) 

G = 200  E2  , 


where 


V 


2 


1 mgal  = 10  5ms  2 , 

(1-3) 

1 E = 1 Ebtvos  = 0.1  mgal/km  = 10  9s~“  . 


This  implies  a curvature  parameter 

X = 7.1  . (1-4) 

The  value  of  Cq  is  taken  from  (Tscherning  and  Rapp, 
1974).  The  correlation  length  £ corresponds  approximately 
to  Kaula's  (1959)  c.f.;  it  would  be  highly  desirable  to  have 
a reliable  new  determination  using  the  presently  available 
data . 

The  value  of  Gq  is  probably  the  weakest  in  the  set 
(1-2).  It  seems  to  be  in  rough  agreement  with  the  scarce 
information  on  gravity  gradients  found  in  various  publication 
(cf.  Schwarz,  1976,  p.15)  and  also  with  the  values  of  gradients 
of  deflections  of  the  vertical  given  by  Bursa  (1974).  Further 
determinations  of  G^  are  urgently  necessary. 

In  their  fundamental  work,  Tscherning  and  Rapp  (1974) 
have  given  a very  complete  discussion  of  possible  analytical 
models  for  global  covariance  functions  and  have  fitted  a model 
to  given  degree  variances  c3  through  cnQ  and  to  the 
variance  C . However,  the  correlation  length  5 comes  out 

O 

too  small  and  G , or  x > i s far  too  large, 
o 

In  the  present  report  we  shall  try  to  fit  a global 
c.f.  model  to  both  the  local  parameters  C , £ , G and  to 

O O 

the  lower  degree  variances.  The  numerical  values  (1-2)  will 
serve  only  for  the  purpose  of  an  example;  the  method  given  is 
applicable  also  to  other  numerical  values. 

The  reason  why  \ for  the  Tscher n i ng - Ra pp  model  is 
too  high  is  that  this  model  is  essentially  a logarithmic  c.f. 
(cf.  Moritz,  1976,  especially  p.48).  Since  the  reciprocal 


3 


distance  c.f.  (ibid.  ,p.35)  has  a \ of  only  3,  we  might  try 
a linear  combination  of  a reciprocal  distance  c.f.  and  a 
logarithmic  c.f.: 


C(P.Q)  = ^ y oJ+2Pn(cos^)  + a 2 l ^+2Pn(cos*)  . (1-5) 
n=3  ‘ n=3 


Here  C ( P , Q)  expresses  the  covariance  between  the  gravity 
anomalies  Ag  at  the  points  P and  Q in  space;  * is  the 
angle  between  the  radius  vectors  r and  r leading  from 

the  earth's  center  of  mass  to  P and  to  Q ; 


r r 

? Q 


and 


r r 

P Q 


(1-6) 


and  Ri  denoting  suitable  constants,  as  are  a and 
a n ; and  the  P^  ( cos-i ) are  the  Legendre  polynomials.  The 
notation  follows  that  of  (Moritz,  1976  , pp. 34-51),  which  may 
also  be  consulted  for  mathematical  details. 

The  disadvantage  of  the  form  (1-5)  rests  in  the  fact 
that  this  model  cannot  be  used  for  the  covariance  function 
K(P,Q)  of  the  anomalous  potential  T because  it  implies 


n+  1 


n + 1 


K(  p ,Q)  = aRji- — P (cos*)  + a,R^ 

a 1 /„i\-  n 


(n-1) 

and  these  series  cannot  be  summed. 

So  we  shall  use  the  form 


-“n(n-l)'  n 


- P (cos*)  , 


(1-7) 


C ( P , Q ) 


n-1 


^?HTrr~Pn(cos*)  + a2^  ( n-7 ) ( n + B ) 


r2P  (cos>) 


(1-8) 


4 


which  has  the  same  local  behavior  since  this  local  behavior 
is  mainly  determined  by  the  higher  degree  terms  for  which 

n - 1 , n- 1 1 

n + A "*  ’ ( n - 2 ) ( n + B ) "*  n 

as  n -*•  * . The  corresponding  K(P,Q)  can  be  summed,  as  well 
as  C(P,Q)  itself,  as  we  shall  see  in  the  next  section. 

2.  The  Model 

Thus  our  covariance  model  will  be 

C ( P , Q ) = a i C t ( P > Q ) + a2C2(P,Q)  (2-1) 

with  ( A and  B always  denote  integers) 

Cjfp.Q)  ■ ’ <2"2> 

C2<P'<»  ’ i(n-Z^6)-jr2p-(c°s'»  ■ <2-3> 

At  sea  level  we  have  in  the  spherical  approximation 
r = r = R , where  R = 6371  km  is  a mean  radius  of  the 

p Q 

earth,  and  a and  a as  defined  by  (1-6)  become  constants. 
Thus  our  model  contains  six  parameters 

“i  ’ a2  ’ 0 i ’ a2  ’ A » 8 (2-4) 


5 


which  may  be  used  to  fix  the  data. 

The  corresponding  c.f.  for  T is  given  by 


K(P,Q)  = a,  Kx  ( P , Q ) + a2K:  (P,Q) 


(2-5) 


with 


= R1  ^ ('n~~i~j  ( n + A )a  1 + 1 Pn  (CQS'f  ) > 


(2-6) 


K.(P.Q)  = Klr^r 

- 3 ^ ' 


n-2) (n+Pic?  1 Pn (cos^ ) 


(2-7) 


If  A 4 -1  , B 4 - 1 and  4 -2  , then  all  these 
series  can  be  summed  in  closed  form  by  a decomposition  into 
partial  fractions,  as  shown  by  Tscherning  and  Rapp  (1974); 
in  fact,  C,(P,Q)  is  their  Model  4 and  C (P,Q)  reduces 
for  A = -2  to  their  Model  3. 

We  have 


n-1  _ , A+l 
n+A  ” n+A 


(2-8) 


n-1  1 , 1 . B + l , 

( n- 2 ) ( n + B ) B + 2 n - 2 n + B , 


(2-9) 


. JL/  i - J-) 

A+ 1 ' n - 1 n + A ' 


n + A ' , 


(2-10) 


n^T  + 
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with 


We  then 


With  Tscherning  and  Rapp  we  put  (ibid.  , p . 3 1 ) 


F (a  »i|» ) = Ion+1P  (t)  , 


(2-12) 


F.(o,<.)  = lian+1Pn(t)  for  i > 0 , 

o 


(2-13) 


%<■>.*)  * 1 f°r  1 S ° 

1 - i 


(2-14) 


t = COSti) 


(2-15) 


have  for  A and  B > 0 


(P.Q) 


■J  F(o1  ,ll>)  - o 1 - CT^t  - =’^P2(t)|  - 


2t  3 
a « a . t 


- (A  + l)ax  Fa(«j1,i|»)  * 7T  ' ATT  - xS?P2(t)'  ' 


(2-16) 


C -,  ( P » Q ) = B+?a  2F-2  ^a2  + 


B+1  Ft  \ 2 

+ B+7"  2 F B ° 2 * ^ ' B~ 


r;t 


b+T  " F+7 


rrp->(t)  ; 


(2-17) 


r: 


! ( p » q ) = a+t  _ a i p 2 ^ 


r: 


°i  ait 


ATT  ' Fa<V*>  ‘ 7T  ' ATT  ' AT7P2(t)  ’ 


(2-18) 


7 


K2  ( P ’ Q)  - g+7F- 2 ^°2  ^ } 


BTT  | F-i  <a2’*)  " °2P2(t)]  + 


TB+r)-(B+r)‘  r 3 ^ J 2 


rFB(o,.*>  - ^ - 


2 

a;  t 


B B + l 


BT?P2 ( t 
(2-19) 


Values  of  A and  B < -3  are  impossible  because 
otherwise  a denominator  in  (2-2),  (2-3),  (2-6)  or  (2-7)  would 
be  zero.  The  values  B = -1  and  -2  are  excluded  because  in 
this  case  the  series  (2-3)  or  (2-7)  cannot  be  summed.  This 
excludes  all  negative  values  for  B ; for  A only  the 
negative  value  -2  is  possible. 

Thus  we  have  for  A = -2  : 


C.(P,Q)  = o l F ( a t ) - al  - Ojt  - a^?2(t):  + 


+ F-2  (J1  ,V>  ) ’ 


(2-20) 


Kx (P,Q)  = - R;  i F_1(a1,b)  - a^P2 (t)  + 


+ rJf_2(0i,*)  • 


(2-21) 


Of  particular  importance  for  the  present  report 
will  be  the  case  A = B = 0 , for  which 


Cx  (P,Q)  = a1  | F ( a x ,*)  - - ojt  - a ^ P 2 ( t )J  - 

- aJ  F o(olf«)  - ojt  - P2  (t)J  , (2-22) 

C 2 ( P » Q ) = 7a2F-2  ^a2  ^ + 

+ \°2  [Fo(02**>  ' °2t  ' ?°2P2(t)  ; {2'23) 

Kx  (P.Q)  = [F_1(a1  ,*)  - 0 i P 2 ( t rj  - 

' R1  [F0(°l»+)  “ ~ fJP2(t>]  ’ (2'24) 

K2  ( P » Q ) = |R2F-2(a2^)  ‘ 

“ R2[F-l(°2**)  * CT2P2(t)]  + 

+ 7R2  [Fo(!V’)  ‘ °2t  ' ?°2P2(t)]  • (2_25) 


L 


Closed  expressions  and  a recursion  formula  for  F 
and  the  F for  i = -2,  -1,  0,  1,  2,  ...  may  be  found  in 

i 

(Tscherning  and  Rapp,  1974).  From  this  report  we  take: 


F ( a , ^ , (2 

F0(or,<|/)  = aln|  , (2 


F_1(a,ti-<)  = a(M  + a t*1  n-^- ) , (2 

F_2(a,*)  = £j(1  + 3ot)M  + c 3 P,  ( t ) 1 n-jj  + ^3(l-t2)  , 

(2 


where 


L=|//l-2ot  + a3  , (2 

M = 1 - L - at  , (2 


N = 1 + L - at 


(2 


For  the  F with  i = 1,  2,  3,  ...  we  have  the 
recursion  formula 


F 


- ( i - 1 ) e”  1 F 


i+1 


f L + ( 2 i - 1 ) tF 


(2 


10 


which  is  valid  for  i > 3 , starting  with 


F 


l 


1 n ( 1 


+ 


2j 

1 -g  + L 


) , 


(2-34) 


F,  = a- 1 ( L - 1 + tFt)  . 


(2-35) 


In  the  same  manner,  closed  expressions  can  also  be 
found  for  the  corresponding  covariance  functions  for  other 
quantities  of  the  anomalous  gravity  field  (Tscherning  and 
Rapp,  1974),  especially  for  second-order  gradients  (Tscherning, 
1976).  We  shall  not  give  them  here,  because  for  the  present 
report  we  need  only  Gq  , for  which  an  approximate  expression 
will  be  derived  in  the  following  section  and  exact  formulas 
in  the  Appendix. 


3.  Variances 

The  variance  Cq  represents  the  value  of  the  co- 
variance  function  C(P,Q)  for  the  case  that  the  points  P 
and  Q coincide  with  a point  at  sea  level,  that  is,  for 


R , 


'p  = 0 . 


11 


From  (2-1)  we  get 


Co  = a 1 C 1 O + a2c20 


(3-1) 


where  C 1 Q and  C0Q  are  the  variances  corresponding  to 
C,(P,Q)  and  C ^ ( P , Q ) , respectively. 

By  substituting  p = 0 in  (2-22)  and  (2-23)  we  get 
for  the  case  A = B = 0 : 


= a1'F(a1>  0 ) ' ' ~ x ' 


1 ! Fo(<v  c)  ' ai  ■ iai  1 ’ 


(3-2) 


' 20  = 2a2F-2(a2’  0)  + 7a2  ! Fo(a2’0)  ' 


2 1 3~j  , , , * 

z 2 ' 2 : • ' 3 ' 3 ' 


The  expressions  (2-30)  through  (2-31)  become  for  p = 0 


Lo  = 5 ’ 


Mo  - 0 • 


(3-4) 


N = 2;  , • 

o 


where  we  have  put 


c = 1 


(3-5) 


Using  these  values  in  (2-26)  through  (2-29)  and 
substituting  in  (3-2)  we  get 


13 


all  ?n(0)  = 1 , we  get  from  (2-12) 


F(J,0)  = la 

o 


n + 1 _ a _ a_ 
1-a  " l 


and 


F(a,0)  - 


a - a 


3 r n 
J - 

4 


30  30 


4 r n 
a ) a 


By  integrating 


h 

i 


n- 1 _ 1 


with  respect  to  a we  get 


V 1 n ,1  .1 

) —a  = 1 n-= - I n— 

“ , n 1-a  <, 

n=  1 


On  changing  the  summation  index  by  putting 


m = n - i 


we  have 


1 = 0 i y = Ini  . 


tn=  1 - i 


m+  1 


i . m + 1 
n=  1 - i 


Therefore  (we  write  again  n in  the  place  of  m ) 


1 n*  1 
-a 


n + 1 


j 1 ~ 1 1 ni 


(3-11) 


(3-12) 


(3-13) 


(3-14) 


(3-15) 


J 


14 


L 


In  this  way,  (2-14)  becomes 


F.(a,0)  = j1  iln—  for  i < 0 
L 4 


whereas  for  (2-13),  with  i > 0 , we  have 


00  CO  — ]_  — 1 

1 = 1-')  (1  = 0 by  definition)  , 

o 1-i  1-i  o 


so  that 


1-i-.  -1 


- 1 


1 n + 1 


F.(o,0)=a  ln—  - 7 r0n  for  i > 0 

i C . “ . n+ i 


1-i 


Thus,  for  A = -2  , we  get  from  (2-20) 


clo  = — + ojln^  , 

l°  ;1  1 M 


wheras  for  A and  B > 0 we  find 


10 


= — - (A+l  )cd~Al  n-1  + ( A + 1 ) 7 it 
'l  1 'l  Cj"* 


2 a 


n+ . 


n + A 


'20 


J2+(B+1)J2~B1  1 B+ 1 
BT2 1n~  ' 


0 n + 2 
2 


4 2 ' + B 


(3-16) 


(3-17) 

(3-18) 

(3-19) 

(3-20) 


For  A = 0 this  reduces  to  (3-6)  and  (3-7),  as  it  should  be. 

In  a similar  way  we  could  readily  compute  the  variances 
K10  and  K , Q for  A ^ 0 and  B/0,  but  we  shall  not 
need  them  in  this  report. 


We  finally  consider  the  gradient  variance  G ^ . 

It  is  the  variance  of  the  mixed  hor i zonta 1 -ver t i ca 1 gradient 
3ag/3x  , or  half  the  variance  of  the  vertical  gradient 
3 A g / ) z ; cf.  (Moritz,  1976).  In  fact,  expanding  C(s)  , the 
c.f . of  yg  , in  powers  of  the  horizontal  distance  s , 


C(s)=a  -as~+as4  -+ 

cl  2 


(3-21) 


and  differentiating: 


C'  (s)  = - 2a  s + 


C " ( s ) = - 2a:  + . 


we  have 


r~  | 

c"(0)  = I |c'(s)j  = - 2aL  , 


so  that  by  eq . (1-31)  (ibid.,  p.10): 


var(jAa) 


- C " ( 0 ) - iic'(s)"1 


- 4ai  - 2Go 


since 


Go  = 2ai 


(3-22) 


by  eq.  (3-9),  ibid.  , p.24. 

Expanding  2.g  in  spatial  spherical  harmonics, 


16 


r , R . n + 2 , 

= 1(7)  *gt 


and  differentiating  along  the  radius  vector  r (with  which 
the  local  z-axis  coincides)  we  get 


3Ag  _ ^n+2/R^n+3 
jr  - i R 'r' 


R V' 


<^9. 


so  that  the  degree  variances  for  3Ag/3r  arise  from  those 
of  Ag  by  multiplication  by  the  factor 


/n+lr 


Hence  we  have  corresponding  to  (2-1) 


G = a G + a.G0 
0 1 10  2 20 


(3-23) 


where  by  (2-2)  and  (2-3) 


J oo 

/-*  l r ^ " 1 / n , 9 \ ^ n+  1 

Gio  = 7r(5n(nt2>  J1 


J ^ 00 

ry 


n - 1 


2 n+1 


G:o  = ^(n-'2)(n+B)(n  + 2)  a2 


(3-24) 

(3-25) 


the  factor  1/2  expresses  the  fact  that  G is  half  the 
variance  of  3Ag/3r  . 

Since  the  gradient  is  a very  local  quantity,  only  the 
hi gher-degree  terms  have  much  influence,  so  that  only  the 
asymptotic  behavior  of  the  terms  in  the  sums  is  essential. 
Hence  we  may  approximately  replace 


J 


17 


L 


i^(nt2)2  by  n(n-l)  , 


n - 1 


rn~2)Tn'  + 6)  ^n  + 2)"  by  n 


which  behave  asymptotically  (for  n - ~)  in  the  same  way. 
We  thus  get,  to  a very  good  approximation, 


J , °o 


1 O 


= —“jl n ( n - 1 ) a" 

on--. 


n + 1 


2fT  3 


(3-26) 


a , oo 


20  2 R 2 3 


n a 


n+  1 


(3-27) 


These  series  can  be  readily  summed.  By  differentiating 


O = 

o 


(3-28) 


twice  with  respect  to  a we  find 


r „ n-  1 
na 


(1-a) 


(3-29) 


n ( n - 1 ) j n 


(1-a) 


(3-30) 


We  thus  have 


r „ n+  1 

nj  , 


. 1 


(1-0 


13 


£ n ( n - 1 ) a : 


n+  1 


2a' 


( 1 : ) 


o • c- 
2 a , = — - 


the  approximations  are  justified  since 


and 


s2 


will  be 


l 


and 


very  sma 11,  so  that 

terms  1/q2  and  2 / z,]  >>  1 

1 


are  very  close  to  1 and  the 


In  this  manner  we  finally  obtain 


l o 


1 


(3-31) 
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2R2C? 


(3-32) 


These  expressions  are  very  simple  but  quite  accurate:  it  is 
not  difficult  ot  verify  that  they  have  a relative  error  of 
order  z,  . Since  c will  be  seen  in  the  next  section  to  be 
around  0.001  to  0.006,  these  expressions  are  accurate  to 
within  a few  percent. 1! 

A thorough  check  of  (3-31)  and  (3-32)  is  provided 
by  expanding  the  c.f.  into  a series  (3-21)  and  using  (3-22). 

In  doing  this,  we  shall  follow  the  procedure  and  the  notations 
of  (Moritz,  1976  , sec. 4). 

By  means  of  (3-5)  and  on  putting 


\ = 2s  in^- 


(3-33) 


= a\ ' 


(3-34) 


the  expressions  (2-30),  (2-31),  and  (2-32)  are  readily  trans- 
formed into 


l) 


Rigorous  expressions  for 
Appendix. 
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and 


20 


will  be  found  in  the 


L = |/c“+p' 


(3-35) 


2 7 1 2 

N = ? + '/z,  + P + 2°  » 


(3-36) 


, 2 2 1 2 

M = c - /c  +o  + 2"d“  . 


(3-37) 


These  expressions  are  still  rigorous. 

We  now  substitute  these  expressions  into  (2-26) 
through  (2-29),  and  these  into  (2-22)  and  (2-23).  We  now  expand 
into  a power  series  with  respect  to  o , neglecting  again 
a relative  error  of  ; . After  somewhat  lengthy  calculation 
we  obtain 


ci(»>  - C10  • * 0(>J)  ' 


(3-38) 


C2(P)  = C 20 


1 2 . n,  4, 

— -o  + 0(c  ) • 


(3-39) 


The  plane  horizontal  distance  s is  related  to  p 


s = 2Rs  i n|-  = R\  = R ( 1 + c ) 


so  that,  with  a relative  error  of 


(3-40) 


This  is  substituted  into  (3-38)  and  (3-39),  so  that 


Cr(s)  » C 


■T — rS  - + . . . 
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C2^S)  C20 


1 s2  + 


4R2c 


whence  (3-31)  and  (3-32)  follow  in  view  of  (3-22) 


4.  Numerical  Model  Fit 


We  shall  now  show  how  the  present  model  can  be  fitted 
to  the  data  by  an  appropriate  choice  of  parameters. 

The  degree  variances  correspondi ng  to  the  model  (2-1) 
through  (2-3)  are  given  by 


c 

n 


n-1  n+2  . n-1  n+2 

H+7T3!  ai  + ( n-2)  ( n+B) °2  a2  • 


(4-1) 


Since  estimates  of  c , through  c _ . (or  even  higher)  are 
available  and  since  each  n from  3 to  20  gives  one  equation 
(4-1),  it  might  be  thought  possible  to  determine  all  six 
parameters  a , a , a , a , A,  B in  this  way.  However,  it 
turns  out  very  soon  that  this  is  impossible  since  the  equation 
system  for  this  purpose  is  very  poorly  conditioned. 

In  fact,  it  is  reasonably  feasible  only  to  determine 
a and  a7  in  this  manner.  The  values  of  , o0,  A,  B must 
be  found  by  other  considerations  and  given  in  advance.  Then 
a and  a,  can  be  obtained  by  a least-squares  fit  to  the 
degree  var i ances  . 

Such  a least-squares  fit  by  means  of  a two  - pa rame ter 
model  has  already  been  done  by  Tscherning  and  Rapp  (1974,  p.23). 
Since  the  "fitted"  degree  variances  are  very  smooth,  it  is 
possible  only  to  take  two  of  them,  say  c,  and  c,  , to 
determine  a 1 and  a0  . Eq.  (4-1)  gives 


A 
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2 5 2 5 

3 +Aa  l 01 1 3 + B"  2a2  “ C 3 • 

(4-2) 

19  _ 2 2 , 19  22 

20+AJ l al  1 8 ( 2 0 + B ) a 2 a 2 ” C 20  • 


If  these  two  equations  are  solved  for  a and  a,  and  the 
values  so  obtained  are  used  in  (4-1)  to  compute  the  inter- 
mediate , it  will  be  seen  that  also  these  c will  agree 

well  with  the  "fitted"  Ts cher n i ng -Ra pp  values. 

From  Tscherning-Rapp  (1974,  p . 23 ) we  thus  take 


c 3 =31.5, 


c 


20 


10.2  . 


(4-3) 


These  are  the  "global"  parameters  to  be  used.  For  the 
parameters  we  take  and  Gq  from  (1-2): 

Cq  = 1800  mgal2  , 

G = 200  E2  ; 

O 


" local" 


(4-4) 


the  correlation  length  will  not  be  considered  at  present. 

To  these  four  constants  (4-3)  and  (4-4)  we  shall  try 
to  fit  the  model  (2-1)  through  (2-3)  with  A = B = 0 , which 
comprises  four  parameters  o^,  o0,  c^,  . These  four 

parameters  will  be  determined  from  the  four  given  values 
c3>  c0o>  Cq,  G,  by  an  iterative  procedure. 

With  A = B = 0 , the  system  (4-2)  becomes 


2 5 

3 a 1 3 1 


+ 


2,5 
3 C 2 a 2 


( 4-5a  ) 


It 


19  22 

70° 1 ai  + 


19  22 


*2  = C20  • 


It  may  be  given  the  simpler  form,  putting  a2/0i  = q * 


ai  + q a2  = fi  ’ 


A 1 _2  2 - 
“l  * T5q  “2  ■ f 2 


where 


TC3  ’ 


22C20 


The  solution  for  a 1 and  a 9 may  be  expressed  as 


a2  = (qi  - y^22)"1^  - f2)  , 


ai  - fi  - q5^ 


Further  equations  are  furnished  by  (3-6)  and  (3-7) 


r - 1 Ju1  1 4 
Clo  ' ~ " Vn“  ' Ti 


. 1 2,u  2,,1  1 3 1 4 

'20  " 7a2^  a2l1n~  ~ 2 " TJ2 
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where 


= 1 - 


= 1 - 


(4-11) 


Then  the  variance  of  the  gravity  anomalies  is  given  by 


C 

O 


a l C l o + 


a2C20 


(4-12) 


For  the  gradient  variance  we  have  (3-31)  and  (3-32): 


G 


10 


1 


G 


20 


(4-13) 


(4-14) 


and 


G 

O 


a 1 G 1 O + *2G20 


(4-15) 


These  are  the  basic  formulas  for  our  iteration,  using 
the  numerical  values  (4-3)  and  (4-4). 

First  we  note  that  ^ and  a,  as  resulting  from 
the  solution  of  ( 4- 5a , b ) are  not  very  sensitive  with  respect 
to  a,  and  . At  any  rate,  a1  and  a.,  will  be  close  to 

(and  slightly  smaller  than)  unity.  We  therefore  start  with 
some  approximate  values  for  a and  a2  , if  they  are  not 
available,  we  may  use  a = a.,  = 1 as  our  starting  values. 
With  j ^ and  a 2 we  get  f,  and  f,  from  (4-7)  and  a . 
and  a from  (4-8). 
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Now  we  can  try  to  determine  better  values  of  o1  and 
a 2 in  such  a way  as  to  get  the  deseired  Cq  = 1800  and 
= 200  . This  may  be  done  in  the  following  manner.  Start 
with  some  a1  , compute 


i.  ■ 1 ■ 'i 


and  G1q  by  (4-13).  Now  c2  can  be  calculated  in  such  a way 


that  the  sum 


aiG10  + a2G20 


gives  the  correct  value  G = 200  . In  fact,  then  G„  must  be 

o 2o 


G20  = “21 (2°°  ‘ aiGlO}  ’ 


(4-16) 


and  5 = 1 - a2  can  be  found  by  solving  (4-14)  for  c 2 • 

With  these  values  o and  a 2 we  can  compute 
by  (4-9),  (4-10)  and  (4-12).  In  general,  the  C thus  ob- 

o 

tained  will  not  be  equal  to  the  desired  value  1300.  We  can, 
however,  repeat  the  computation  of  C for  various  assumed  a 

O 

always  finding  the  corresponding  a by  (4-16).  If  in  this 

way  we  get  values  for  C that  are  partly  < 1800  and  partly  > 

o 

1800,  we  may  interpolate  the  argument  o for  the  desired 
value  C = 1800  . Again  the  corresponding  a.  can  be  com- 

o 2 

puted  using  (4-16)  . 

Thus  both  G and  C are  fitted,  but  the  values 

o o 

a and  a2  are  not  yet  completely  correct.  We  thus  repeat 
the  computation  of  a and  on  the  basis  of  the  values 

o and  a n last  obtained.  Using  these  values  ^ and  a 0 
we  compute  better  a and  cn  by  again  assuming  a suitable 
a (usually  it  is  best  to  take  the  last  value),  computing 


L 
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the  correspond  i ng  a2  by  (4-16)  and  evaluating  the  corresponding 
C . For  a slightly  different  assumed  a the  process  is 
repeated,  and  the  corresponding  to  the  desired  C_  = 1800 

i s i nterpol ated  . 

Now  we  can  anew  compute  and  a2  and  repeat  the 

iteration  as  long  as  necessary.  The  principle  is  always  the 
same:  assume  » compute  the  corresponding  a2  in  such  a 

way  as  to  ensure  the  correct  G = 200  , calculate  C for 

O O 

different  and  interpolate  a correspondi ng  to  the 

desired  Cq  = 1800  . The  iteration  is  necessary  because  a 
and  a 2 depend  on  g and  a 2 , but  fortunately  this  depen- 
dence is  rather  slight,  so  that,  in  general,  the  iteration  will 
converge  well. 

For  the  given  values  (4-3)  and  (4-4)  the  result  is 

j1  = 0.993  963  , (4-17) 

g 2 = 0.999  281  ; 

ax  = 9.907  , (4-18) 

a 2 = 37.7  7 . 

The  correlation  length  £ of  the  c.f.  with  these 
parameters  can  be  computed  from  the  functional  representa ti on 
(2-22)  and  (2-23),  using  eqs.  (2-26)  through  (2-32).  It  is 
the  horizontal  distance  5 for  which 

C(5)  = |c  = 900  mga 1 2 . (4-19) 


The  computation  may  again  be  done  by  interpolat 
C ( -;,  ) and  C(£2)  for  two  assumed  arguments  ■;  L 
and  interpolate  Z for  which  C(;)  = 900  . If 


ion:  calculate 

and  z 2 > 

■ . and  £ , 
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have  not  yet  been  sufficiently  close  to  5 , the  process  may 
be  repeated.  The  result  is 


C = 68.91  km  , (4-20) 

correspond i ng  to  an  angular  distance 


v = 0.6197°  . 


(4-21) 


The  variance  of  the  anomalous  potential  T is 

obtained  from  (3-9)  and  (3-10): 


K = a,  K + cyK  = 186.3  + 381.8 

o 1 10  2 20 


or 


Kq  = 568.1  ( k g a 1 • m ) ' 


(4-22) 


The  correspondi ng  variance  of  geoidal  heights  N is 


var(N)  = = 591.5  m2  , 

9 


(4-23) 


g = 0.98  kgal  being  a mean  value  of  gravity 


5.  Conclusions 


The  numerical  model  discussed  in  the  preceding  section 

does  not  fit  all  parameters  (1-2);  it  only  fits  C and  G 

00 

Its  correlation  length  ^ has  been  found  to  be  68.9  km,  which 
is  smaller  than  the  value  80  km  corresponding  to  (1-2). 


It  is  possible  to  get  a different  value  of  - by 
selcting  A or  B in  a different  way?  Rather  conclusively 
the  answer  appears  to  be  NO. 

To  see  this,  let  us  fix  the  value 


and 


A * -2 


vary  B . Then  the  system  (4-2)  becomes 


(5-1 


a?ai  + 3TB°2a2  = d 3 ’ 

22  A 1 22  , 

°l  ai  20  + Ba2  a 2 d20  ’ 


(5-2 


where  we  have 

put 

d 3 - 

2C3  5 

d 2 0 

. lie 
19  20  ' 

As  a 

first. 

very 

crude,  approximation  we  put 

3i  = a2  = 1 » 

obtaining 

as  a sol uti on 

(5-3 


*1  4 d 20  - T7(Vd20)(3  + B) 
a2  = "IT(d  3_d20)  (3  + B)  (20  + B) 


(5-4 


With  the  values  (4-3)  this  becomes 


= 8.6  - 0.36  B , 

= 0.36(3+B) (20+B)  . 


(5-5 
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VJ hat  are  the  possible  limits  between  which  B can 
vary?  A 1 ower  limit  is 

B = 0 , (5-6) 

as  we  have  seen  on  p.7.  An  upper  limit  may  be  obtained  by  the 
condition 

<*!  > 0 ; (5-7) 

for  negative  a1  the  positive  definiteness  of  the  covariance 
function  C(P,Q)  can  no  longer  be  guaranteed.  The  limiting 
case  = 0 corresponds,  by  (5-5),  to 

B = 8.6/0.36  = 24  . (5-8) 

This  corresponds  to  the  case  of  Tscherni ng-Rapp 
(1974,  p.22),  in  which  a model  with  only  the  function  C0(P,Q) 
is  used.  However,  for  this  model  G >>200  E“  . This  indicates 
that  if  Gq  = 200  E“  is  imposed,  the  upper  possible  limit 
should  be  smaller.  It  is  found  to  be 

B = 10  ; (5-9) 

for  larger  B the  conditions  C = 1800  and  G = 200  turn 
out  to  be  incompatible. 

In  fact,  for  each  B we  can  apply  the  iterative 
procedure  described  in  the  preceding  section,  to  find  1 1 , 
a,  and  , a2  corresponding  to  the  values  Cq  and  G_ 

given  in  advance.  It  is  found  that  this  is  possible  for  all 
B from  0 to  10,  whereas  for  B = 11  and  higher,  the  values 
of  corresponding  to  G ^ = 200  are  all  < 1300  , so  that 

an  interpolation  of  for  = 1300  can  no  longer  be 

perf ormed . 
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For  B = 0 we  thus  find 

a = 0.994  032  , 

c,  = 0.999  484  ; 

and  the  correlation  length 

£ = 63.7  km  . 

For  B = 10  we  obtain 

a = 0.995  455  , a = 6. 199  , 

j2  = 0.997  953  ; a0  = 127.28  ; 

£ = 60.7  km  . 


= 10.000  , 
= 18.18  ; 


(5-11) 

(5-12) 

(5-13) 


(The  comparison  between  these  values  of  a and  with 

the  approximate  values  given  by  (5-5)  shows  that  the  latter  is 
only  a rough  approximation.) 

It  is  thus  seen  that,  for  A = -2  and  B variable, 
the  possible  values  of  £ lie  between  rather  narrow  bounds: 


60.7  < £ < 68.7  km  , 


(5-14) 


the  value  £ = 30  km  cannot  be  reached  by  this  model  (it 
can  be  expected  that  for  B = 1,  2,  ...,  9 the  corresponding 
values  £ lie  indeed  within  the  limits  (5-14)). 

It  is  seen  that  even  by  selecting  a 
different  A , the  situation  will  not  change  essentially.  The 
model  retains  its  general  character,  the  coefficients  of  C 
decreasing  asymptotically  as 
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(5-15) 


corresponding  to  the  reciprocal  distance  covariance  function. 


and  those  of  C.,  as 


1 n + 2 

na2 


(5-16) 


correspond i ng  to  the  logarithmic  covariance  function.  What  may 
change  are  the  numerical  values  for  the  lower  and  upper  bounds 
in  (5-14);  in  fact,  for  A = B = 0 we  had,  by  (4-20) 

5 = 63.9  km  . 

It  can  be  expected,  however,  that  these  changes  are  small.  At 
any  rate,  we  shall  again  have 


5i  S 5 S 


(5-17) 


l and  C2  differening  perhaps  by  10  km. 

It  thus  appears  that  given  values  of  the  lower  degree 
variances  (say,  from  3 to  20),  together  with  C and  G 

O O 

rather  narrowly  determine  the  correlation  distance  s . For 
instance,  choosing  in  the  present  case 


7('i+52) 


(5-18) 


taking  A = -2  and  fitting  B to  this  i,  , the  correlation 
distance  of  the  model  will  be  wrong  at  most  by,  say,  5 km. 

The  exact  upper  and  lower  bounds  t,  and  5 if 
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both  A and  B vary  may  probably  best  be  found  by  a computer 
experiment  since,  on  the  one  hand,  the  formulas  for  general 
A and  B are  not  difficult  to  program  but,  on  the  other  hand, 
the  required  hand  calculations  are  very  laborious.  If  this  is 
done,  then  the  exact  expressions  for  Gio  and  G?q  given  in 
the  Appendix  might  be  employed  rather  than  the  approximate 
ones  used  above. 

A condition  of  form  (5-14)  seems  to  be  quite  generally 
true,  not  only  for  the  special  model  (1-8)  used  in  this  report. 
Even  the  "general"  5 and  should  not  differ  very  much 

from  those  for  the  model  (1-8)  with  variable  A and  B , and 
even  from  those  for  A = -2  and  B variable. 

For  this  we  may  adduce  the  following  argument.  Other 
possible  covariance  models  have  coefficients 


n a 


n 


2 n + 2 
a 


„ 3 n+ : 
n a 


etc . , 


but  these  are  all  ruled  out  because  they  would  have  increasing, 
rather  than  decreasing,  degree  variances. 

The  simplest  model  with  decreasing  degree  variances 
is  (5-16),  corres pond i ng  to  a logarithmic  covariance  function, 
but  such  a model  alone  will  give  too  high  a G ^ ; therefore 
the  reciprocal  distance  model  (5-15)  is  also  necessary.  More 
general  would  be  an  expression  with  degree  variances 


n+2 
a . 


+ 


+ 


+ 


(5-19) 


but  it  can  be  expected  that  the  first  two  terms,  c or  res  pond i ng 
to  the  model  used  in  the  present  report,  will  be  dominant. 

The  expression  (5-19)  represents  a linear  combination 
of  degree  variances 


1 
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n“kan+2  with  k = 0,  1,  2,  ...  (5-20) 

Even  non-integer  values,  say  k = 0.5  , do  not  seem  to  have 
an  essentially  different  behavior. 

Such  a different  behavior  could  only  be  expected  if, 
in  an  essential  way,  the  degree  variances  behave  capriciously, 
for  instance  decreasing  from  n = 3 to  30  , then  increasing 
to  n = 100  , etc.  If  they  decrease  monotoni cal ly , at  least  on 
the  average,  then  the  behavior  can  be  expected  to  be  similar 
to  our  present  covariance  model. 

This  fact,  if  it  is  true,  is  of  very  grea  t i mportance . 
Already  Cq  , G^  and  the  lower  degree  variances  together 
determine  to  correlation  length  5 rather  well.  Some  values, 
such  as  1 = 80  km  with  the  presently  assumed  numerical  data, 
are  excluded.  This  emphasizes  the  importance  of  determining 
a good  value  for  the  gradient  variance  G . This  quantity, 
which  is  a very  local  one,  is  the  same  for  local  and  global 
covariance  functions;  this  is  not  true  for  t,  . This  G , 

O 

rather  than  ^ , provides  a basis  for  calculating  a global  c.f. 
from  local  covariance  functions. 
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APPENDIX 

Gradient  Variances 

Let  us  calculate  the  difference  between  the  exact 
expressions  for  the  gradient  variances  G and  G2Q  , given 
by  (3-24)  and  (3-25),  and  the  approximate  formulas  (3-26)  and 
(3-27). 

We  use  the  following  auxiliary  formulas  which  can  be 
directly  verified: 


£^(n  + 2)2  - n(n-l)  = (4-A)n  + A(A-3)  - (A  + ^~2): 


n-2  n + E 


-( n + 2 ) - n = 5 - B + 


i + 2 ) n-i 


(B  + l) 

( B_  2 ) 

(n  + B) 

( A- 1 ) 


The  difference  between  (3-24)  and  (3-26)  is  therefore 


.exact  .approximate 

G10  ' Gio 


2 


jl  ! S^(n  + 2)2  - n ( n - 1 )~|  a?+1  . 


2R2 3 


~2 

j (4-A)[na”+1  + A(A-3 
2R  L 3 1 3 


- ( A+l ) ( A-2 ) 


n + 1 “1 

°1  I 

n + A 


(A-2) 
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In  view  of  (3-28)  and  (3-29)  we  have 


I na  i 


n+  1 


(1-0,) 


2 9 3 

2 ' °1  ' 2ai  ’ 


( A - 3 ) 


r n+ 1 4 / •>  2 , 

lal  = ox (1  + a1  + o:  + . . . ) 


l‘al  ’ 


( A - 4 ) 


and  by  (3-15), 


r 1 n+l  1 - A -1  1 

- °i  lnI^ 


1 - A 


so  that 


r 1 n+l 
3n+Aai 


a 


1 -A 
1 


1 


n+  1 
ai 

n+7T  5 


(A-5) 


for  A = -2  , the  value  of  the  last  sum  is  taken  to  be  zero. 
Thus  ( A - 2 ) becomes,  with  1 - a = c 2 , 


2 

(4-A)(4  - a*  - 2aJ)  + 

1 1 

1 

4 

+ A(A-3)-^  - 

9 an+1  "1 

- (A  + l)(A-2)2(aJ-Al  nl—  - j -ij- ) i . ( A- 6 ) 

The  difference  between  (3-25)  and  (3-27)  is  transformed 
similarly.  We  have 


A G , 


10 


2R" 


gexact  _ ^approximate 
20  " 20 


2R23  L 
2 


(n-1) 

1 1 

ln  + £( 

Li 

( n-2) 

i ( 

! n + 8) 

i 

n+1  , 1 6 r a 2 


^ / r-  n \ r _ n + l 

_ + BT23-—  + 
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By  (A-4)  and  (A- 5)  this  becomes,  with  1 - a 2 = , 
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The  exact  sums  of  the  approximate  expressions  (3-26) 
and  (3-27),  are,  by  the  last  formula  on  p .17  and  the  first  on 
P-  18, 


^ap proximate  1 "2  9 3 

>|  • 

(A-3) 

2 r 2 

rt  i 

r apor ox i 
G20 

mate  _ J 2 2 2 3 

2R2  ; 2 2 2 

'-2 

( A - 9 ) 

We  thus  have  the 

exact  expressions 
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the  summands  on  the  right-hand  side  being  given  by  (A-8)  and 
(A-6),  and  by  (A-9)  and  (A-7),  respectively. 

The  corrections  aG1q  and  are  order  ; , 

so  that  the  formulas  (3-31)  and  (3-32)  are,  in  fact,  accurate 
to  that  order. 
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